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; Abstract 

■ Global existence results in the past time direction of cosmological mod- 

■») | els with collisionless matter and a massless scalar field are presented. It 

is shown that the singularity is crushing and that the Kretschmann scalar 
diverges uniformly as the singularity is approached. In the case without 
Vlasov matter, the singularity is velocity dominated and the generalized 
V-^. Kasner exponents converge at each spatial point as the singularity is ap- 

proached. 



X ' 1 Introduction 



In ^21 the mathematical study of inhomogeneous cosmological solutions of the 
Einstein equations coupled to the Vlasov equation and a linear scalar field was 
begun. A local existence theorem and continuation criteria were proved for 
solutions with surface symmetry (i.e., spherical, plane or hyperbolic symmetry). 
In ^1] this was used to prove a global existence theorem in the expanding 
direction. In the case of plane symmetry with only a scalar field it was shown 
that the solutions are future geodesically complete. In the present paper we 
turn to the study of the past time direction, i.e., the approach to the initial 
singularity. 

In the case where the matter is given by the Vlasov equation alone certain 
results on the initial singularity have been obtained by Rein JJJ . He showed that 
if the maximum momentum of the particles remains bounded on any interval 
where the solution exists then, for spherical or plane symmetry, the solutions 
can be extended up to t = 0. In the case of hyperbolic symmetry he showed 
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that the corresponding result holds provided the metric function fi (defined be- 
low) is initially negative. Weaver jTB] showed that the solution always exists up 
to t = in the case of plane symmetry. This is a special case of a result she 
proved for solutions of the Einstein- Vlasov system with T 2 -symmetry. Tchap- 
nda |12j extended the result from plane symmetry to spherical symmetry and 
to hyperbolic symmetry under the assumption that fj, is initially negative or the 
assumption that fi is bounded on any interval where the solution exists. His 
results also allow for a negative cosmological constant. 

In the case where the matter is given by a scalar field alone existence up 
to t = for solutions with plane symmetry follows from the work of |1U| . In 
that case it is only necessary to solve a linear hyperbolic equation which is 
identical to the essential field equation in polarized Gowdy spacetimes. Using 
results of [5] it was shown that the initial singularity is a curvature singularity 
(the Kretschmann scalar blows up uniformly there) and a crushing singular- 
ity (the mean curvature of the hypersurfaces of constant t blows up uniformly 
as t —> 0). For matter described by the Vlasov equation alone the results al- 
ready mentioned can be combined with theorems in [7] to show that t = is a 
crushing singularity where the Kretschmann scalar blows up uniformly. For an 
open set of initial data belonging to a restricted class the asymptotic behaviour 
can be described more precisely. For convenience these will be referred to as 
small data. The generalized Kasner exponents converge uniformly to the values 
(2/3,2/3,-1/3). This means that in a certain sense the solution is approxi- 
mated near the singularity by a Kasner solution with those particular values of 
the Kasner exponents. 

In this paper some of the above results will be generalized to the case where 
both Vlasov matter and a scalar field are present. It is shown that the solutions 
exist up to t = for spherical and plane symmetry and in the case of hyperbolic 
symmetry when fj, satisfies the restrictions mentioned above. The Kretschmann 
scalar and the mean curvature of the hypersurfaces of constant t blow up uni- 
formly as t — > 0. Under the assumption that the maximum momentum of any 
particle in the radial direction decays like a positive power of t as t —> it is 
possible to obtain further interesting estimates. Unfortunately we did not suc- 
ceed in generalizing the small data results for the Vlasov equation alone to the 
present case. For a scalar field alone these estimates allow the detailed asymp- 
totics to be determined for all solutions. The generalized Kasner exponents 
converge to the values ((1 — o(r))/2, (1 — a(r))/2, a(r)) for a continuous function 
a(r). It may be noted that these solutions have the same type of singularity as 
those constructed from data on the singularity in pp . 

Let us recall the formulation of the Einstein- Vlasov-scalar field system as 
shown in and |15|. We consider a four-dimensional spacetime manifold M, 
with local coordinates (x a ) = (t, x 1 ) on which x° — t denotes the time and (x 1 ) 
the space coordinates. Greek indices always run from to 3, and Latin ones 
from 1 to 3. On M, a Lorentzian metric g is given with signature (— , +, +, +). 
We consider a self-gravitating collisionless gas and restrict ourselves to the case 
where all particles have the same rest mass, normalized to 1, and move forward 
in time. We denote by (p a ) the momenta of the particles. The conservation 



2 



of the quantity g a pp a p requires that the phase space of the particle is the 
seven-dimensional submanifold 

PM = {g af 3P a P p = -1; p° > 0} 

of TM which is coordinatized by (t,x l ,p l ). If the coordinates are such that 
the components goi vanish then the component p° is expressed by the other 
coordinates via 

P° = \f Z 9^°\l 1 + flyPV 

The distribution function of the particles is a non-negative real- valued function 
denoted by /, that is defined on PM. In addition we consider a massless scalar 
field 4> which is a real- valued function on M. The Einstein- Vlasov-scalar field 
system now reads: 

V Q V Q = 

f i dp dp dp 1 
T a p = - I fp a p/3 \ g 1 1 1- (V Q </>V/3</> - -g a pV u4>V v 4>) 

JR 3 PO * 

where p a — g a /3P^ , \g\ denotes the modulus of the determinant of the metric 
9a0, the Christoffel symbols, G a p the Einstein tensor, and T a p the energy- 
momentum tensor. 

Note that since the contribution of / to the energy-momentum tensor is 
divergence-free , the form of the contribution of the scalar field to the energy- 
momentum tensor determines the field equation for (j>. 

We refer to (HJ for the notion of spherical, plane and hyperbolic symmetry. 
We now consider a solution of the Einstein- Vlasov-scalar field system where all 
unknowns are invariant under one of these symmetries. We write the system in 
areal coordinates, i.e., coordinates are chosen such that R = t, where R is the 
area radius function on a surface of symmetry. The circumstances under which 
coordinates of this type exist are discussed in 0. In such coordinates the metric 
g takes the form 

ds 2 = _ e Mt,r) dt 2 + e 2Ht,r) dr 2 + f 2 {dQ 2 + ^2 (1 ^ 

where 

{s'm8 for k = 1 (spherical symmetry); 
1 for k — (plane symmetry); 

sinh# for k = — 1 (hyperbolic symmetry) 

t > denotes a time-like coordinate, r £ R and (8, if) range in the domains 
[0, 7r] x [0, 2tt], [0, 2tt] x [0, 2tt], [0, oo[x [0, 2tt] respectively, and stand for angular 
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coordinates. The functions A and [i are periodic in r with period 1. It has been 
shown in that due to the symmetry, / can be written as a function of 



t, r, w := eV and F := t 4 [(p 2 ) 2 + sin 2 k 0(p 3 ) 2 } 



i.e., / = f(t,r,w,F). In these variables, we have p° = e~^\J\ + w 2 + F/t 2 . 
The scalar field is a function of t and r which is periodic in r with period 1. 

We denote by a dot and by a prime the derivatives of the metric components 
and of the scalar field with respect to t and r respectively. Using the results of 
the complete Einstein- Vlasov-scalar field system can be written as follows: 



dtf+ 7= ^—— d r f-(\w + e> 1 - x »'y/l + w*+F/t*)d w f = (1.2) 
\J 1 + ur + b J t 2 

e- 2 ^(2t\ + l) + k = 8nt 2 p (1.3) 

e- 2 ^(2t(i -l)-k= STrt 2 p (1.4) 

fi' = -4:TTte x+ Vj (1.5) 

e ~ 2 V + M '( M ' - A')) - e- 2 »(\ + (A + -j-)(A - fi)) = 4nq (1.6) 

e- 2 V - e" 2 ^ - e- 2 "(A - /t + -)</» - e~ 2A (A' - //)<£' = (1.7) 
where <|1.7|l is the wave equation in <f> and : 

7T /' + 00 /" + 00 

p(t, r) = e- 2 "Too(t, r) = ^ J J Vl + w 2 + F/t 2 f(t, r, to, F)dFdw 

+ i(e^ 2 + e - 2 V 2 ) 

(1.8) 



* J-oo Jo \J\ + w 2 + F/t 2 



+ I( e -2^2 + e -2A^2 ) 



(1.9) 



. +oo /-+oo 

j(t,r) = -e- {x+ ^T 01 (t,r) = - / / wf(t,r,w,F)dFdw - e -< A+ ">M' 

— oo JO 

(1.10) 
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t 4 J-oo Jo v/l + ^ + F/i 2 ' V ' 

(1.11) 

We prescribe initial data at time t = 1: 

f(l,r,w,F) = f(r,w,F), A(l,r) = A(r), M (l, r) = £(r), 

o 

0(1, r) = 0(r), 0(1, r) = ■0( r ) 

The choice t = 1 is made only for convenience. Analogous results hold in the 
case of prescribed data on any hypersurface t = t > 0. 

The paper is organized as follows. In section 2, we show that the solution of 
the Cauchy problem corresponding to system l|1.2fl - l|l.llfl exists for all t G]0, 1]. 
In section 3, we analyze the asymptotic behaviour of solutions as t — > 0. The 
paper ends with a discussion of some interesting open problems. 

2 Global existence in the past 

We use the continuation criterion in the following local existence result. 

Theorem 2.1 Let f e C^M 2 x [0,oo[) with 

f(r + 1, w, F) = f(r, w, F) for (r, w, F) E K 2 x [0, oo[, / > 0, and 

o 

w := sup{|iw||(r, w, F) e supp/} < oo 

o 

Fo := sup{i 7 '|(r, w, F) e supp/} < oo 
Let X,ip e C^R), C 2 (R) with A(r) = A(r + 1), £(r) = /i(r + 1), 

o o 

0(r) = <j>(r + 1), — i ) { r + 1) a7l <^ $1.5\) satisfied for t = 1. JTien </iere 

exists a unique, left maximal, regular solution (/, A,/x, 0) of system ^1.^ - ^1. 11^1 

° ° o ° 

(/, A, /u, 0)(1) = (/, A,/z, 0) and 0(1) = tp on a time interval ]T,1[ with 

Te[0,l[. // 

i. sup{|t/;||(t, r, w, F) € supp/} < oo , 

£ sup{(e~ 2 ^0 2 + e" 2A 0' 2 )(t, r); r6K}< oo, 

5. /i is bounded, 

then T = 0. If k > or /j, < then condition 3 is automatically satisfied. 

This is the content of theorems 4.4 and 4.5 in ^3]. For a regular solution, all 
derivatives which appear in the system exist and are continuous by definition 
(see H3|). 

In order to obtain the global existence of solutions, we prove the following results 
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Lemma 2.2 Let D+ = e~^dt + er x d r ; D~ = e~^d t - e~ x d r ; 



X = <j>e-» - 4>'e~ x 



b = 



b'e- x 



,'„-A 



. V / c= (-A- i)e^ + M 'e 

and define X 2 — e^X , Y 2 = e M Y. Then as a consequence of the field equations 
X and Y satisfy the system 



D+X = aX + bY 



(2.1) 



D~Y = bX + cY (2.2) 
If in addition the field equations M.ty - fT^]) are satisfied, then X 2 and Y 2 satisfy 



D+X 2 = e"[- - *irt(j> - p)]X 2 - —Y 2 



D-Y 2 = -^X 2 + e^* - Avt{p - p)]Y 2 



Proof: This results from a straightforward calculation. □ 
Lemma 2.3 Define X 2 and Y 2 as in lem,ma Wl\ and let 
B(t) = sup{(|X 2 | 2 + |Y 2 | 2 ) 1/2 (^) ; rel} 
l(t) = S up{j + e 2 ^+4Trt(p-p)}(t,r) ; r 6 R} 

If (X 2 ,Y 2 ) is a solution of f2.S\l - ^4]) , then we obtain the estimate 
B{tf < B(l) 2 + 2 f l(s)B(s) 2 ds 



with t e ]T,1], T > 0. 

Proof: We deduce from system I|2.3 p - I12.4[l : 

~k 



D + Xl = 2e» 



t 



Ant(p — p) 



Xi - 2—X 2 Y 2 



D~Yo 



-2—X 2 Y 2 + 2e^ 



4irt(p — p) 



Y 2 



(2.3) 
(2.4) 



(2.5) 



On the corresponding characteristic curves (t, 7j), i = 1, 2 of the wave equation, 
(see |Tjl]) D + or D~ is equal to e"^^ and then 



|xf(t, 7l (t)) = 2e 2 " 



k 

- ~ 4nt(p - p) 



Xl(t, ll (t))--X 2 Y 2 (f ll (t)) 



dt 



Y 2 2 {t,l2{t)) = --X 2 Y 2 (t, l2 (t)) + 2e 2 » 



4:Trt(p — p) 



G 



Integrate each of the two previous equations on [t, 1] and obtain respectively : 
X|(t,7 1 (t))=Jf|(l, 7l (l)) + 2 f{e^[--+ATTs{p-p)]Xl + -X 2 Y 2 }{ Slll { S ))d S 

.It s s 



<X|(l, 7l (l)) + 2^\[l + e ^(-^+47r S (p-p))]X| + lr 2 2 }( S , 7l ( S ))d S 



r/(t, 72 (i))<^(i, 72 (i))+2 / {^2+ 



1 i 



1 k 

— + e 2 "( +47rs(/9-p)) 

2s s 



Y 2 2 }(s, l2 (s))ds 



Add the two previous inequalities and take the supremum over space to obtain 
estimate (|2.5|) .D 

Unless otherwise specified in what follows constants denoted by C will be pos- 
itive, may depend on the initial data and may change their value from line to 
line. 

Proposition 2.4 Let (f,\,/j,,<j>) be a solution of the full system on 
a left maximal interval of existence ]T, X], T > 0, with initial data as in theorem 

H3 // 

1. Q(*)=sup{M|(r,ic > /)esiipp/(t), te]T,l]}<oo 

2. fj, is bounded. 

then T = 0. If k > or /j, < then condition 2 is automatically satisfied. 

Proof: We need to prove that K(t) = sup{(|X| 2 + \Y\ 2 ) 1 / 2 (t, r) ; r e M} is 
bounded for all i G]T, 1]; where X and Y are defined in lemma l2~2l Subtract 
the two equations (|1.8|l - H1.9fl to obtain : 

p-n= — / / — fdldw 
t'J-ooJo ^/l + w 2 +F/t^ 



/ + F/t 2 fdFdw 

-oo JO 

C 

< -g ; since < oo 

Using (4.15) of which shows that e 2fl = O(i), we obtain : 

,/ \ c |fe| 

Then (|2.5|) implies : 

B(t) 2 < B{1) 2 + 2 f (C(l + -) + -)B(s) 2 ds 
Jt s s 



(2.6) 
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And by Gronwall's lemma, 

B(t) 2 < cr c - 2 

We have from l|1.9|) , 

/-+0O r + oo 2 1 

P(s,r)<- / ^ /(s , r , w ,F)dFd«; + -e^(X 2 2 -(-y 2 2 )( S ,r) 

S J-oc JO M 4 

< 5 + ;k 2Mi? ^ 2 

Then using p.4|l . we obtain the estimate : 



i 8 ^ f 2 / v . 
fc / s p(s, rjds 



<Cr 1 {\ + s- c e- 2ll ds). 



By Gronwall's lemma, we deduce that e 2 '' < Ci 1 exp[y^p(l — i 1 c )]. 
Therefore, 

(X 2 + Y 2 )(t, r) = e- 2 ^(X 2 2 + F 2 2 )(i, r) < Cr 3 ~ c exp[ I ^(l - i 1 " 6 ")] 

i.e K(t) is bounded. We conclude by theorem 12. II that T = O.D 

In the next theorem, we prove that the solution exists on the full interval ]0, 1]. 



Theorem 2.5 Consider a solution of the Einstein- Vlasov- scalar field system 
with fc > and initial data given for t = 1. TTien ttis solution exists on the 
whole interval ]0, 1]. //fc < cmd // < 0, £/ie same result holds. 

Proof : The strategy of the proof is the following : suppose we have a solution 
on an interval ]T, 1] with T > 0. We want to show that the solution can be 
extended to the past. By consideration of the maximal interval of existence this 
will prove the assertion. 

Firstly let us prove that under the hypotheses of the theorem, /i is bounded 
above. From the field equation (|1.4|l and since p(s, r) > 0, we have for fc > 0, 

t 

For the case fc = -1, e- 2f * > e [- 1 + 1 > 1 which gives the upper bound of /z 

o 

for /i < 0. In either case, /i is bounded above. 
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Now, let us prove that w is bounded. Consider the following rescaled version 
of w, called Ui, which has been inspired by the works of ^B] (p. 1090) and 
(p. 5336): 

Ui = —w. 
2t 

If we prove that [i is bounded below then the boundedness of u\ will imply the 
boundedness of w. So let us show that p is bounded below under the assumption 
that u\ is bounded. We have 

^(te- 2 ») = -k-8irt 2 p. (2.7) 
at 

Transforming the integral term defining p to u\ as an integration variable instead 
of w yields 

/ = Ul fdFd Ul + -(X 2 +Y 2 ); 

where X and Y are defined in lemma l2~21 The integrand in the first term in p 
can then be estimated by 47re~ 2 ' 1 |wi|. We have from ll.8fl - H1.9fl : 

m, s * f°° f°° (l + F/t 2 )2td l ,, , 
e 2 »{p-p) = -, / .= K = o l =d= ==fdFdu l 



t 2 J-oo Jo y/1 + 4i 2 e" 2 ^ 2 + F/t 2 ' 
<— / {l + F/t 2 )e^fdFd Ul 

where u\ is the maximum modulus of u\ on the support of / at a given time. We 
can then estimate from (12 .5|) l(s) by h(s) = Csup{l+s -1 +s~ 2 e M '0i(s, r); r E M} 
and B{t) 2 by 5(l) 2 exp( // h(s)ds). Thus 



X 2 + Y 2 < e- 2fl B(t) 2 < e- 2fl B(l) 2 exp (J 



h(s)ds 



Therefore, using the bound for \i and u\, p can be estimated by Ce 2ai and so 
(|2.7|l implies that 

\j t (te^)\<C(l + te- 2 n, 

Integrating this with respect to t over [i, 1] and using the Gronwall inequality 
implies that te~ 2 ^ is bounded on ]T, 1]; that is fj, is bounded below on the given 
time interval. 

The next step is to prove that u\ is bounded. To this end, it suffices to get 
a suitable integral inequality for u\. Since u\ — u\(t, r(t)), we can compute iii : 

e M g^ 
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i.e. 



u\= ( p + rp - - j ui + —w (2.8) 



We have 

u! = -Airte^j, f = 

^/l + w 2 + F/t 2 

and 

. l _|_ fc g 2M 

w = 4nte 2 ' i {jy / l + w 2 + F/t 2 - pw) H — w 

so that 12. implies the following : 



u'i|ui = e 2 ' 1 



-4nt(p -„) + £] Ul \ Ul \ + 2^^=£±£gP_ (2.9) 
t\ y/\ + 4t 2 e- 2 ^u 2 + F/t 2 



In order to estimate the modulus of the first term on the right hand side of 
equation l|2.9|l . we need an estimate of e 2/i (p — p)u 2 . For convenience let log + 
be defined by log + (a;) = log a; when logcc is positive and log + (a;) = 0, otherwise. 
Then estimating the integral defining p — p shows that 

p-p<C(l + log + (w)), 



i.e. 



P -P < C(l + log + (ui) - p). 



The expression —p, is not under control ; however the expression we wish to 
estimate contains a factor e 2 ^. The function p i— » —pe 2 ^ has an absolute maxi- 
mum at —1/2 which is (l/2)e _1 . Thus the first term on the right hand side of 
equation (|2.9|) can be estimated by Cu\{l + log + (rii)). 

Next the second term on the right hand side of equation (|2.9|) will be esti- 
mated. By definition 

7 ' ' - ■ - ■ ■ ,/„ n A _ 



3 = p J J o wf(t, r, w, F)dFdw - ^' e~^~ x = j x + j 2 

The first term of j can be estimated by Cw 2 , i.e. 

\]i\<Cu 2 e- 2 » 

and the second term < -^e~ 2,1 B(t) 2 ; so that it suffices to estimate the 

quantity 

(,l + B {tm + F/t 2 ) 
y/l + 4t 2 e- 2 "u? + F/t 2 

in order to estimate the second term on the right hand side of equation 12. 9|) . 
But since p and i -1 are bounded on the interval being considered, the quantity 
H2.10JI can be estimated by C(u 2 + B(t) 2 ). Thus adding the estimates for the 
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first and second terms on the right hand side of (|2.9[) allows us to deduce from 
that 

MM < Cu\(l + log + (tii)) + C{u\ + B(t) 2 ) 



i.e 

||ki| 2 |<C(S 1 ) 2 (l + log + (u?)) + C J B(<) 2 
Integrating over [t, 1] gives : 

u\{t) < vll(l) + C J [u\{s){\ + log + (S 2 ( S ))) + B(s) 2 ] ds (2.11) 

We deduce from the estimate of p — p and from inequality l|2.5[l . that l(s) can 
be estimated by C(l + log + (ui)). We then obtain 



B(ty < B(iy + C j (1 + log + (^)(s))B(s)^s (2.12) 
Adding (|2.11(l and Ij2.12jl gives estimate : 

u\{t)+B{tf < u\{l)+B{lf+C (l+ul+B(s) 2 ) [1 + log + (l + ul(s) + B(s) 2 )] ds 
Set v(t) = u\{t) + B(t) 2 , then the above estimate can be written 

v(t) < v(l) + C [ (1 + «(*)) [1 + log + (l + «(«))] ds (2.13) 



By the comparison principle for solutions of integral equations, it is enough to 
show that the solution of the integral equation 

a(t) =a{l)+C J (l + a(s))(l + log + (l + a(s)))ds 

is bounded. The solution a(t) is a non-increasing function. Thus either a(t) < e 
everywhere, in which case the desired conclusion is immediate or there is some 
T\ in ]T, 1] such that e < a(t) on ]T, T{\ . We take T\ maximal with that property. 
Then it follows on ]T, T{\ that inequality 

a(t)<C[l+ [ a(s)(l + log o(s))eLs] 



holds for a constant C. The boundedness of a(t) follows from that of the solution 
of the differential equation x = Cx(l + logx) which is exp(exp(Ci) — 1). In 
either case a(t) is bounded. Thus u 2 and B{t) 2 are bounded i.e w and K (t) are 
bounded. The proof of the theorem is complete using proposition 12 .41 □ 
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3 On past asymptotic behaviour 

In this section we examine the behaviour of solutions as t — > 0. Firstly we 
generalize the work of Ringstrom S310-S311) to bound the quantity 

\cj)'\e IJ, ~ x by C|ilogt| _1 , where C is a positive constant. 

Lemma 3.1 Let A x = \{-(j> + j^- t + 4> 1 e^ x ) 2 and 

A 2 = |(— </> + tl Q gt — tj>'e^~ x ) 2 with t e]0, 1[. If (f> satisfies the wave equation, 
then 

(d t + e^d r )A x = -1(1 + + -f-f + 0' 2 e 2 ^ 2A ] 

At log t t log t 

+ ^(1 + ^-)0' 2 e 2 ^ 2A + 1(A - A + -)(* - <^- A )(-0 + tt^7 + <^~ A ) 

2i logi 4 t ilogt 

(3.1) 

(ft + e ^ dr)M _ _ _L (1 + _L )M + + ^»-«] 

+ 5" + 1^"^ + j (a - * + > + ^ )( -* + Tsa - ^ 

(3.2) 



Proof : This results from a straightforward calculation. □ 



Proposition 3.2 Let (/, A, /1, 4>) be a left maximal solution of the Einstein- 
Vlasov- scalar field system on the interval ]T, 1], < T < e _1 . Assume that 

Q(t) = sup{| w I |(r,u>,F) e supp/(t)} < Ci Q 

/or some positive constants C, a and for some t G]T, e -1 ]. TTien 

(-0+-^) 2 + 0' 2 e 2 ^ 2A < C(t logi)- 2 (3.3) 
rlogr 

Proof : Consider the two characteristic curves (t, 7i(i)) and (t, 72 (i)) of the 
wave operator. Since t e]0, e _1 ], the term (1 + i^j)4>' 2e ' 2fJ '^ 2X is nonnegative 
and (-0 + j^- t ) 2 + <j>' 2 e 2 ^- 2X = 4(Ai + A 2 ), then from (fSHJ : 



(dt + e^dM^t^t)) > -1(1 + 7^7)^1 + A 2 )(t, 71 W) 

- J(A - A + l)(-0 + 0'e^ A )(-0 + 77^7 + <^- A )(*, 7i W) 
Similarly, we deduce from l|3.2l) that : 

(9 t - e^ x d r )A 2 (t, l2 {t)) > -1(1 + T-l^i + 4,)(f, -*(*)) 

- 1(A - A + 1)("^ ~ 0'^ A )(-0 + 77^7 - 0'^ A )(t, 72 W) 



(3.4) 



(3.5) 
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Since Q(t) < Ct a , we can bound p - p by Ct" 3+Q (see e 2 ^ 1 < Ct; then 

subtracting (|1.3I) - (11.41) gives 

(A + i = -^-+47rte 2 "(/5-p) < C{l+t- 1+a ); and from (£3, 1(a) can 
be bounded by + C + Cs~ 1+Q . We deduce from i|2.5|) (consider the integral 
term in the interval [t, e -1 ]) that 

B(t) 2 < B^e- 1 ) 2 exp[2 Jf 1 (s- 1 + C + Cs- 1+a )ds] i.e., B(t) 2 < Ct~ 2 . 
Therefore \<j>(t)\ and \(f)'\e^~ x (t) are bounded each by Ci _1 . We can then have 
a lower bound of the second term of the right hand side of each inequality (|3.4|) 
and which is -C(t- 2 + t- 3+a ). Then 

(d t + e^dM^t^t)) > -hi + r^-MA, + A 2 )(t,^(t)) - C(t- 2 + t- 3+a ) 

t logr 

and 

(fit - e^ x d r )A 2 (t, 72 (i)) > ~\(l + ^)( A i + A 2)(i, 72 W) - C(t- 2 + t- 3+a ) 

On the corresponding characteristic, we have dt + d r or dt — e' 1 fir equal 
to Jr. Take the supremum in the space of each of the above two inequalities 
and add them. Then 

±(A X + A 2 )(t, r) > -- (1 + J-)(^ a + A 2 )(i, r) - C(i- 2 + ^ 3+ ") 
at r log r 

Set u(t) = (Ax + A 2 )(t) and v(t) = (t\ogt) 2 u(t). If v(t) is bounded, then we 
conclude that u(t) is bounded by C(t log £)~ 2 . Let us prove that v(t) is bounded. 
We have : 

dv , , . 9 , du 2 , 1 . . 

^- (ilog<)( ^ + i (1 + ioi7 )M) 

> -C(log<) 2 (l + ^ 1+Q ) 

Then, u(t) < i^e" 1 ) + C ff\l + S - 1+a )(log s) 2 ds < v^' 1 ) + C. We obtain 
the desired conclusion of the proposition. □ 

In the case / = 0, we obtain from the previous proposition, theorems 12.51 
and 12. II the global existence of solutions and the above estimates hold on the 
whole interval ]0, e -1 ]. In general we do not know how to use this proposition 
to obtain precise asymptotics. It seems that if the estimates could be improved 
slightly they would allow a bootstrap argument on the bound for Q similar to 
that used in |7j . These estimates have been included here in the hope that they 
might help someone else to complete the argument. 

Next we prove that the curvature invariant R a pjsR al3lS called the Kretschman 
scalar blows up as t — ► 0. Then there is a spacetime singularity and the space- 
time cannot be extended further. 

Theorem 3.3 Let (/, A, /i, (f>) be a regular solution of the surface- symmetric 
Einstein- Vlasov- scalar field system on the interval ]0, 1] with data given for 
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t = 1. Then 

{R al3lS R aMS ){t,r) > - (infe^ + fc) 2 , (3.6) 

Proof We can use the following expression for the Kretschman scalar from [7]. 
R aPl sR aplS = 4[e- 2 V + m'(m' - A')) - e-^i'X + A(A - /i))] 2 

+ l(e-^ + fe) 2 



The first term K\ is nonnegative and can be dropped. Inserting the expres- 
sions 



e-^A = kirbp - 



k + e~ 2fl 
2t 



e 2m /i = 4-7rtp 



fc + e~ 2 ^ 
2T 



; e 



- A ~V = -47rfj 



into the formula for K2 yields 



16-K 2 r{p z +p 2 - 2j z ) - 4nt{p-p) 



k + e' 2 ^ (k + e- 2 ^) 2 



2t 2 



Now 



00 />oo 



< 



/>oo />oo />oo />oo 2 

-S[ / ^l + W 2 + F/t 2 fdFdw]H / -__==/dFd«;]* 
1 J -00 Jo J -00 Jo A/l + tf z + F/t /i 



b 2 e 2/i + 0' 2 e 2A ) by the Cauchy-Schwarz inequality. 



In fact the above inequality holds in general for all choices of matter satisfying 
the dominant energy condition. Therefore 



K 2 > 



> 



t 2 

4 
J 2 



8irt 2 (p~p) 2 -Ant(p-p) 



k + e- 2t " {k + er 2 ^f 



2t 2 



4att(p — p) 



k + e 



> 0. 
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Recalling the expression for e 2m we get 



~ 2 ^ + fc = - — + —f s 2 p(s,r)ds 



e~ 2fl + k inf e~ 2/i + k 



and thus 



K 3 = -(e-^ + k) 2 >-(irfe-^ + k) 
We obtain (|3.bj) and deduce that 

\im(R al37S R a ^ s )(t,r) = oo, 

uniformly in r £ M. 

Next we prove that the singularity at t — is a crushing singularity i.e. 
the mean curvature of the surfaces of constant t blows up. In the case where 
there is only a scalar field and no Vlasov contribution this singularity is velocity 
dominated i.e the generalized Kasner exponents have limits as t — > 0. 

Theorem 3.4 Let (/, A, fi, <f>) be a regular solution of the surface- symmetric 
Einstein- Vlasov- scalar field system on the interval ]0, 1] with initial data given 
on t = 1. Let 

K{t,r) := -e - " (\(t,r) + j 
denote the mean curvature of the hyper surf aces of constant t. Then 

K{t,r) < -Cr 3 / 2 , 
where C is a positive constant. 

Proof We use the same argument as in [7] and obtain the following : 
K(t, r) = -(A + 2) e -/* ; A = e 2 ^ (4wtp - ^f^) > -e 2 ^ ( h± ^ L ) and 

k-3e- 2 » „ 



K(t,r) < 



2t 



For fc = or fc = — 1, 



and the estimate 



K(t,r)<-- ( 



-2m > + k 



implies 

K (t,r) < _A ( infe ' 2M + fc ) i/2 < _ Ct -V2 where C= -(infe-^ + fc) 1 / 2 . 
2£ t 2 
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For k = 1 we have 



thus 



-2„ > £ > 1 = fc 



P 2M _ Q 

tf(t,r) < (— g— ) — 

o 

e _A1 infe _/i 
< < 



t ~ t 3 / 2 
<-Ct- 3/2 where C = infe^. 



We deduce from above that 



uniformly in r G 



lim K(t, r) = —oo, 
t->o v y 



Theorem 3.5 Let (A, /x, 4>) be a regular solution of the Einstein- scalar field sys- 
tem with spherical, plane or hyperbolic symmetry on the interval]®, 1] with initial 
data given at t = 1 . Then 

, K}(t,r) , . , K$(t,r) , iffft.r) 1, 

lim V - / = ofr) ; lim 2 , ' / = lim 3 ; - / = -(1 - air)), 

uniformly in r £R, where 

K\{t,r) K 2 {t,r) K 3 (t,r) 
K(t,r) ' if(t,r) ' #(t,r) 

are i/ie generalized Kasner exponents and a(r) a continuous function of r. 

Proof We have as in [7] 

gjfcr) = t\(t,r) K 2 (t,r) = K 3 {t,r) = 1 
K{t,r) tk(t,r)+2' K(t,r) K{t,r) t\{t,r)+2 

As we have seen previously 

e 2 M (t,r) < ct 

which implies that 

e 2M*.r) _> o as t-> 
Let t G^e" 1 ] and i e]0,i o ]. From lET3ll . 

log i log t t log I 
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Since t 1 (log t) 2 is integrable on the interval (0, to) it follows that we can define 

A(r) = lim ^l = - / \]ags)-\fo,r) - far)/ slog s)ds 

t^o logt log to Jo 

Since from E3, (> - ■ n ^ i ) = 0((t| logtl)" 1 ), we have 

i^-^ + Ofdlogtl)- 1 ) 
logt 

so that t(j> — > A(r) as £ — s- 0. Inequality 1)3. 3[) shows also that 
(f)' 2 e 2fl ~ 2X = 0((t\ logt|)~ 2 ). Using these limits, we have 

t\(t,r) = 2n(t 2 <j} 2 +t 2 cj)' 2 e 2 ^- 2X ) - -e 2 » - ~ -> 2ttA(V) 2 - i as t -> 0, uniformly in r. 
We take a(r) = I^TrpT^ *° com pl ete the proof. □ 

4 Discussion and outlook 

It is an open problem to remove the restriction on /i in Theorem 2.5, even when 
the scalar field vanishes. The only example where existence up to t = is 
known to fail is a vacuum solution, the pseudo-Schwarzschild solution (cf. the 
discussion in jHJ, p. 115). Perhaps it can only fail in the vacuum case. As 
mentioned in the introduction, the result for the plane symmetric case extends 
to solutions of the Einstein- Vlasov system with T 2 -symmetry |16| . 

Once existence up to t — is known the ideal goal is to obtain detailed 
information on the asymptotics. For the Einstein- Vlasov system this has been 
done in [7] for small data. In Proposition 3.2 an analogue of some parts of the 
proof of the small data theorem in .7. is obtained but this is not sufficient in order 
to determine the asymptotic behaviour. There is a formal similarity between 
the Einstein- Vlasov-scalar field system with plane symmetry and the Einstein- 
Vlasov system with polarized Gowdy symmetry. If the problem of asymptotics 
could be solved for the first problem it would probably lead to valuable insights 
for the second problem. 

In the case of general Gowdy symmetry, results on asymptotics are available 
in the vacuum case but they are hard to obtain JT]. In the yet more general 
case of T 2 -symmetry the only thing known is a construction of a class of vacuum 
solutions with prescribed asymptotics g|. 

Acknowledgements : The authors acknowledge support by a research 
grant from the VolkswagenStiftung, Federal Republic of Germany. The ma- 
jor part of this work was completed during a visit of D. Tegankong to the 
Max-Planck Institute for Gravitational Physics, Golm, Germany. 

References 

[1] L. Andersson and A. D. Rcndall 2001 Quiescent cosmological singularities. 
Commun. Math. Phys. 218, 479-511. 



17 



[2] H. Andreasson, G. Rein and A. D. Rcndall 2003 On the Einstein- Vlasov 
system with hyperbolic symmetry. Math. Proc. Camb. Phil. Soc. 124, 529- 
549. 

[3] J. Ehlers 1973 Survey of General Relativity Theory. In: W. Israel (ed.), 
Relativity, Astrophysics and Cosmology, 1-125, Reidel, Dordrecht. 

[4] J. Isenberg and S. Kichenassamy 1999 Asymptotic behaviour in polarized 
T 2 -symmetric vacuum spacetimes. J. Math. Phys. 40, 340-352. 

[5] J. Isenberg and V. Moncrief 1990 Asymptotic behaviour of the gravitational 
field and the nature of the initial singularity in Gowdy spacetimes. Ann. 
Phys. (NY) 199, 84-122. 

[6] V. Moncrief and D. Eardley 1981 The global existence problem and cosmic 
censorship in general relativity. Gen. Rel. Grav. 13, 887-892. 

[7] G. Rein 1996 Cosmological solutions of the Vlasov-Einstein system with 
spherical, plane, and hyperbolic symmetry Math. Proc. Camb. Phil. Soc. 
119, 739-762. 

[8] G. Rein 1995 The Vlasov-Einstein system with surface symmetry Habilita- 
tion thesis, Munich. 

[9] A. D. Rendall 1995 Crushing singularities in spacetimes with spherical, plane 
and hyperbolic symmetry. Class. Quantum Grav. 12, 1517-1533. 

[10] A. D. Rcndall 1995 On the nature of singularities in plane symmetric scalar 
field cosmologies. Gen. Rel. Grav. 27, 213-221. 

[11] H. Ringstrom 2004 Asymptotic expansions close to the singularity in 
Gowdy spacetimes Class. Quantum Grav. 21, S305-S322. 

[12] S.B. Tchapnda 2004 Structure of solutions near the initial singularity for 
the surface-symmetric Einstein- Vlasov system Class. Quantum Grav. 21, 
5333-5346. 

[13] D. Tegankong, N. Noutchegueme and A.D. Rcndall 2004 Local existence 
and continuation criteria for solutions of the Einstein- Vlasov-scalar field sys- 
tem with surface symmetry Journ. Hyperb. Diff. Eq. Vol 1, No. 4, 691-724. 

[14] D. Tegankong 2005 Global existence and future asymptotic behaviour for 
solutions of the Einstein- Vlasov-scalar field system with surface symmetry 
Class. Quantum Grav. 22, 2381-2391. 

[15] D. Tegankong 2005 Cosmological solutions of the Einstein- Vlasov-scalar 
field system PhD thesis Technische Universitat Berlin 

[16] M. Weaver 2004 On the area of the symmetry orbits in T 2 symmetric 
spacetimes with Vlasov matter Class. Quantum Grav. 21, 1079-1098. 



18 



